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We consider the Hopf-Oleinik normal derivative lemma for elliptic and parabolic equations under 
minimal restrictions on lower-order coefficients. Boundary gradient estimates for solutions are also 



established. 

1 Introduction 

Qualitative theory of partial differential equations is in intensive development over last half of 
century. In this paper we discuss the Hopf-Oleinik Lemma, one of the most important tools in 
studying solutions to elliptic and parabolic equations, in particular, the key argument in the 
proof of uniqueness theorems. 



! For the Laplace operator this property is well known for one hundred years, starting from 

a pioneer paper of Zaremba [Z], and reads as follows. Let dft G C 2 and let C = —A. Then, if 
^ ; G df2, we have 

du 

Cu = f > in Q; u(x) > u(0) in Q —(0) < 0. (ZHO) 

For general operators of non- divergence type with bounded measurable coefficients this 
result was established in elliptic case independently by E. Hopf [HoJ and O.A. Oleinik [O] and 
in parabolic case by L. Nirenberg [M]. Later the efforts of many mathematicians were aimed 
at the reduction of the boundary smoothnesa3- They established that the sharp condition for 
(ZHO) to fulfil is the Dini condition for the boundary normal, see, e.g., [HiJ. In a weakened 
form (the existence a boundary point x° in any neighborhood of the origin and a direction 
£ such that fj(x°) < 0) this fact holds true for a much wider class of domains including all 
Lipschitz ones, see |Naj for elliptic equations and [K] for parabolic ones. Note that all these 
results are related to classical solutions, i.e. u G C 2 (i?) in elliptic case and u G C 2,1 (Q) in 
parabolic case. 

Now let us consider generalized {strong) solutions for non-divergence type equations 

Cu = -Qij^DiDjU + bi(x)D lU = f(x); (NDE) 

Mu = d t u - a,ij(x; t)DiDjU + h(x; t)Diu = f(x; t), (NDP) 

i.e. we assume D(Du) G L n ^ oc (fi) in (NDE) and d t u, D(Du) G L n+1 j oc (Q) in (NDP) (in the 
parabolic case also some anisotropic spaces are admissible). 



*Partially supported by RFBR grant 09-01-00729 and by grant NSh.4210.2010.1. 
^ee also an earlier paper [G] for equations with Holder continuous leading coefficients. 



We always suppose that operators under consideration are uniformly elliptic (parabolic), 
i.e. for all values of arguments 



m 2 < < ^r, ^« n . (i) 

where v is a positive constant. Note that we can also assume = dji without loss of generality. 

The properties of generalized solutions to the equations (NDE)-(NDP), under assumption 
that leading coefficients are only measurable, were investigated in a number of papers^]. The 
problem of our interest is how "bad" may be lower-order coefficients b{ to ensure the Hopf- 
Oleinik Lemma to hold true. We provide sharp conditions for this. We also touch the topic 
closely related to (ZHO), especially in idea of proof, namely, the gradient estimates at the 
boundary. 

Note that for divergence type equations 

-D i (a ij (x)D j u) + bi(x)DiU = 0; (DE) 

d t u — Di(aij(x; t)Dju) + b,i(x\ t)DiU = 0. (DP) 

(ZHO) does not hold. The simplest counterexample is the function u(x) = x\ + 2x2|xi| which 
is positive in the upper half-plane, satisfies the equation (DE) with 

bi = 0; f = 0, 



{Ct-ijj 



1 — signal, 
-sign(xi) 2 



but u(0, 0) = and D 2 u(0, 0) = 0. 

Moreover, even continuity of does not improve the situation. Let us describe correspond- 
ing counterexample^]. 

Let Q be a convex domain, and let G dfi. Assume that at the neighborhood of the origin 
dQ is the graph of a function x n = <f)(x'). Finally, suppose that <fi G C 1 but D'tfi is not Dini 
continuous at the origin. 

As it was mentioned the Hopf-Oleinik lemma for the Laplacian fails in such domain. Now we 
rectify the boundary near the origin and obtain an operator of the form (DE) with continuous 
leading coefficients and 6j = for which (ZHO) fails in smooth domain. Considering functions 
depending only on spatial variables we see that this example works also for the parabolic 
operator (DP). 

The paper is organized as follows. In Section 2 we deal with elliptic equations, Section 
3 is devoted to parabolic equations. In both sections we use the "composite" variant of the 
A.D. Aleksandrov maximum estimate ( |Li00j ; see also [AN95] for a weaker version) and slightly 
modify classical techniques due to Ladyzhenskaya-Ural'tseva |LU88j . see also |S10j . 

Let us recall some notation, x = (x\, . . . , x n -x,x n ) = (x', x n ) is a vector in M. n , n > 2, with 
the Euclidean norm \x\; (x;t) is a point in M n+1 . 

i? is a domain in MJ 1 and <9i? is its boundary; n = (nj(x)) is the unit vector of the outward 
normal to dQ at the point x. 

For a cylinder Q = f2x]0,T[ we denote by d"Q = <9i?x]0,T[ its lateral surface and by 
d'Q = d"Q U {£2 x {0}} its parabolic boundary. 

2 We mention in this connection a quite recent paper |A-Z| discussing some degenerate elliptic equations. 
3 See also [GTj Problem 3.9]. 



We define 



B r (x°) = {x G W 1 : \x - x°\ < r}, B r = B r (0); 

B r ,h(x ') = {xe R n : \x' - x°'\ <r,0<x n < h}; B r>h = B T , h (0) 

Q r (x°- t°) = B r (x°)x]t° - r 2 ; t% Q r = Q r (0; 0) 

Q r , h (x '; t°) = B r , h (x°>) x]t° - r 2 ; t% Q r , h = Q r>h (0 



0). 



The indices i,j vary from 1 to n. Repeated indices indicate summation. 

The symbol Z)j denotes the operator of differentiation with respect to x^, in particular, 
Du = (D\u, . . . , D n u) = (D'u, D n u) is the gradient of u. d t u stands for the derivative of u with 
respect to t. 

We denote by || ■ the norm in L q (f2). We introduce two scales of anisotropic spaces: 
L q/ (Q) = L t (]0,T[^ L q {Q)) with the norm \\f\\ qAQ = ||||/(-;*)|| q ,n|' 



10,T[ 



L q/ (Q) = L q (p -> L e (]0,T[)) with the norm ||/||~ AQ = .)||/,]o,r[| 



q,Q' 



Obviously, L q q (Q) = L q q (Q) = L q (Q). Further, by the Minkowskii inequality, 
\\f\\; M <\\f\\ 9 M for q>£; \\f\\ qAQ < \\f\\~ M for q < 
We denote by L q< e(Q) the space 



L q/ (Q) n L q/ (Q) 



with the norm \\f\\ qA Q = max{||/||^ )Q , \\f\\ qAQ }. 



L q ,e(Q), q> 
L q ,g{Q), q < 



Remark 1. Note that we always deal with the space L q ^(Q) i.e. take the more strong 
of two norms. The reason is that up to now anisotropic versions of the Aleksandrov-Krylov 
maximum principle (see |N87j . |N01j ) are proved only in terms of stronger norm. 

We set /_)_ = max{/, 0}, /_ = max{— /, 0}. 

Following [LiOO], we say that u : [0, 1] — > 1R+ belongs to the class T>\ if uj(1) — 1, u is 
continuous and increasing while oj(a)/a is summable and decreasing. In this case we define 

We use letters M, N, C (with or without indices) to denote various constants. To indicate 
that, say, N depends on some parameters, we list them in the parentheses: N(. . .). 



2 Elliptic case 

Recall that in this section we assume D(Du) G L n ; oc (i7). 

The next statement is a particular case of [LiOOj Theorem 3.2]. 

Proposition 2.1. Let £ be an operator of the form fNDEj in a bounded, strictly Lipschitz 
domain Q, and let the condition (QJ) be satisfied. Suppose also that the vector function h = (6j) 
can be written as follows: 

b = b (1) + b (2) ; |b (1) | eL n {n); (2) 

1-91, ™ w(rf/diam(i?)) ^ 
b (2) <<B—-^ — - v JJ , ujeVi, 
d 

where d = d(x) = dist(x, dfi). 



Then for any solution of Cu = f in Q satisfying u\ dQ < 0, the following estimate holds: 

diam(J?) 

U < N \\f+\\n,{u>0}, 

V 

'provided ||bW|| ni r2 < 03o, where N and 03o depend only on n, v, 05 and the Lipschitz constant 

of on. 

Now we prove a quantitative version of tne maximum principle, the so-called "boundary 
growth lemma" (its versions for |b| G L n are proved, e.g., in |LU85t Lemma 2.5'] and |S10t 
Lemma 2.6]). 

Lemma 2.2. Let C be an operator of the form (~N~D~E) in B p , p , p < R, and let the conditions 
(QP, (TJP and 

| b (2) (x)| <^^nAR), ueVlt (3 ) 

•En 

be satisfied. Suppose also that ||b^|| n R < 03o where 05o = 05o(?i., v i 05) is the constant from 
Proposition 2.1. If u is a nonnegative solution of Cu = f > in B p , p satisfying u > k on 
dB PiP H {x n = 0} for some k > then for £ < ^ the inequality 



u 



> k ■ (P(n, v, 03, u, - N^n, v , 03, u, f) • (||b« \\ n>Bp>p + <Bu(p/R))) (4) 



holds in with some positive constants (3 and Ni. 

Proof. Consider the barrier function 

2 U'|2 



w(x) = (l-A^) + 2(1 + A) (<p(x n /p) - tp{\/A)) - 1^-, (5) 
where (cf. [LlOO] ) 



V?(s) 



(exp (— X^ajj - ljrfcr, 



while A > — — is a constant to be defined later. 

— V 

Direct calculation shows that for x G B n p. 

-OijDiDjW < \ ■ [ - vA 2 - v{\ + A)<p"{r) + (n- l)^ 1 ] < ~ (1 + A)<f?'(r)\ 

^ 1} A^< |b (1) | •-• [i(l-ir) + (l + iy(r) + l] < |b«|--(l+A)(l + ^(r)); 
p 7 J p 

6- 2) Aw < ^ • - At) + (1 + AV (r) + 1] • 05 ^ < (1 + A) ^ (1 + ^( T )) 

p T p T 

(here r = x n /i?). Since ip"(r) = f ^ (1 + y/(r)), we have £u> < d(A, v,<B)\b^\p- 1 . 

Further, w(x) < for \x'\ = p, < x n < 4 and for < p, x n = 4. Finally, w(a;) < 1 for 
'I < P, = 0. This gives few — u < on dB p ^. 



\x 



Proposition 2.1 gives for x E B 



u(x)>kw(x)-kC 2 p-\\(Cw) + \\ n:B > kw(x)-kC 3 \\(bW\\ n , B 



where C 2 depends only on n, v and 33 while C3 depends on the same quantities and on A. 

Now we observe that (p(s) = o(s) as s — > 0. Thus, we can choose A = A(n, u, ( B, cu, £) > ^ n ~ 1 
so large that 2(1 + < |. Since (1 — |) 2 — | — (1 — £) 2 = f (1 - ?) > f , this gives 

'5£ 



u>A;.(-i-C4(7i,i/,?8,a;,0l|b (1) mB p , /> ) + = A; 1 in S (w)p x. (6) 
Now we consider the set K„ = B a a \ B e» . Note that coefficients fej 2 ' 1 are bounded on this 

p P'JA 1 

set, and 

1 

||b« \\ n>Kp < C 5 (n)QS ( I * ^ ^(n, 1/, ®, w, 0®w(p/i2). (7) 



X 

4.4 



We apply "the ink-spot expansion lemma" ( |LU85t Lemma 2.2]) and obtain 

u>k x - [n{n, v, 93, -C 7 (n,u, ||b (1) + b (2) ||„ iis - p j in S (1 _ e)p>(1 _£ )p \ B (W)P) || . 

By © and © we arrive at (@j). □ 

Lemma 2.2'. Lei £ fre as in Lemma 2.2. If u is a nonnegative solution of Cu — f > 
zn $ pp satisfying u > k on dB P)P fl {x n = p} for some k > i/ien /or £ < | i/ie inequality 
/ioWs m ^5(i-e)p,p \ 

Proof. This statement is more simple than Lemma 2.2. Consider the set K p = £> PiP \B p 

Since coefficients fe^ are bounded on this set and ||b( 2 )|| n ^ is under control, we can apply 
standard boundary growth lemma, and the statement follows. □ 

Remark 2. If we replace the assumption / > by /_ G L n (B PtP ), the estimate (jlj) holds 
true with additional term —N 2 (n, u, 23, u, £)p- ||/_ ||n,B p ,p i n the right-hand side. The proof runs 
without changes. 

Now we prove the main result of this Section. 

Theorem 2.3. Let C be an operator of the form fNDEj in Br^, and let the conditions 
(OP, (E)) and be satisfied. Suppose also that for p < R 

\\b^\\ nt B P , P <^Mp/R). (8) 

Then 

1. Any solution of Cu = f < in Br } r such that u\ Xn=0 < and u(0) = satisfies 

sup < — — • sup u, 

0<x n <R/2 X n R B R/2 ,R/2 



Consequently, if D n u(0) exists then (£)„«) + (0) is finite. 



2. Any positive solution of Cu — f > in Br r such that u(0) = satisfies inf u (°' x «) > 

o- /.' 2 Xn 

0. Consequently, if D n u(0) exists, it is positive. If, in addition, f = 0, the following estimate 
holds: 

, nf u(0,x n ) > N _ u(0,R/2) 

0<x n <R/2 x n ~ 3 R/2 



The constants depend on n, v, 03, 03 1, tu and the the moduli of continuity of |b| in 

Proof. We introduce the sequence of cylinders B Pk; h k , k > 0, where p k = 2~ k po, h k = ( k p k , 
while po < R and the sequence ( k I will be chosen later. 
Denote by Mjjr , k > 1, the quantities 

M+ = sup r , i > sup M~ = inf 



Note that in the case 2 M^T > 0. 
We define two function sequences 

i + , <^ + {x n /R) 2 <f-(x n /R) 

v k = u- Kqh k — ; w fc = M fc fr fe _ - u, 

(p + {h k /R) ip {h k /R) 

where, similarly to Lemma 2.2, 

s 

^(s) = J exp (t— Zu{o-))da, (9) 
o 

and denote 

V k = v\\ M k = .] I,. : $ = (p + in the case 1; 

V k = v\] M k = M' k : $ = tp~ in the case 2. 

It is easy to see that V k \ _ n < while the definition of M k gives V k < on the top of the 
cylinder B PkM . 



To estimate V k} we refine a trick from [S10] . Let x° G B Pk -h k ,h k - Assume first that x° n < - 



k 

■ - ' | — | c 3 r' k '"rev "re — 2 

Then we apply Lemma 2. 43 to the (positive) function M k h k — V k in Bh k ,h k (x 0/ ) (with regard to 
Remark 2). This gives for x G Bh^ ^(x ') 

2 ' 2 

M k h k - V k (x) > M k h k ■ [p{n, v, 23, u, 1/2)- 

- iVi(n, i/, 03, w, 1/2) ■ (||b« ||„, Bpfc , Pfc + - 

- jV 2 (n,i/, 03, 1/2)71* • ||(£^) + ||„, Bhfcihfc( xOo- (10) 

We suppose that po/-R is so small that the quantity in the square brackets is greater that |. 
Further, direct calculation similar to Lemma 2.2 shows that the assumptions of theorem imply 

CV^M.lb^'ixjR)^^ in B Pk)hk . 

Note that (p + is concave, <fi~ is convex, and both of them are increasing. Therefore, 

&(x /R) Hk,R < max{1 ' $/(1)} 

9 [Xn/K) Hh/R) ~ $(i) • 

Substituting these inequalities into (ITU]) we arrive at 



V k (x) <M k h k - 



1-0/2 + C 8 {n, v, 03, Lo)\\bW\\ n>B (x o l for ^(x 0/ ). 



lr To proceed we suppose that po/-R is so small that ||b|| nj jB < %$q. 



In particular, this estimate is valid for x = x°. If x° n > we get the same estimate using 
Lemma 2.2' instead of Lemma 2.2. 

Taking supremum w.r.t x°, we obtain 



sup V k < M k h k ■ (1 - (3/2 + C 8 <BMp k /R)) . 

B Pk -h k ,h k 

Repeating previous arguments provides for m < 

sup V k < M k h k ■ ((1 - 0/2) m + C& ^W ) ■ 

Setting m = |_^~J 5 we arrive at 

M k h k / / p k+ i\ u(p k /R)\ 

where A = - ln(l - 0/2) > 0. 
Therefore, for x G B Pk+lt h k 

] k{X \ t < M fc7fc , (11) 

where 7fc = ^ ^ • (exp ( - + C7 8 »x^). 

Estimate (ITTj) implies in the cases 1 and 2, respectively, 

M^^M^iT,), (12) 

where 



, h k ip + (x n /R) h k (p + (h k+l /R) 

q 1 — . Slip — — * ' 

(f+(h k /R) < Xn < hk max{x n ,h k+1 } cp+(h k /R) h k+1 
~_ _ h k _ <p-(x n /R) _ h k _ ip-(h k+1 /R) 

(f~(h k /R) h k+1 <x n <h k x n (f~(h k /R) 



b k+l 



Since lim = 1, we have 

ik=^s 1 



$(V#> ^ $(i)' 



(13) 



Thus, f lT2|) gives in the cases 1 and 2, respectively, 

«£-. s IS • n (1 ± p s S • n (i ^ • mM { i '*( i )})- 

We set Cfc = and choose k Q so large and po/R so small that 71 • max{l, $(1)} < |. Note 
that fc and po/-R satisfying all the conditions imposed depend only on n, u, u, 23 and QSi. 
Now we observe that the first terms in 7^ form a convergent series. The same is true for 



the second terms, since 



00 

00 „ 

^u(2- k po/R) x / uj(2- s po/R)ds-l u (p /R). 
k=i { 



Thus, the series ^ fc 7fc converges. Therefore, the infinite products II ± = rj(l±7fc-max{l, $(1)}) 

k 

also converge, and we obtain in the cases 1 and 2, respectively, 



n ± M 1 ± 



Thus, all M£ are bounded in the case 1, and all M k are separated from zero in the case 2. 
Further, we note that M± < -r- sup u. This completes the proof of the statement 1. 

If / = 0, then we set K = Err \ Br^/2- Similarly to Lemma 2.2, Hb^ + b( 2 ^|| ra ^ is 
bounded. Therefore, by the Harnack inequality ( |S10t Theorem 3.3]), Mf x ■ This 

completes the proof of the statement 2. □ 

Remark 3. If we replace in the case 1 the assumption / < by / = /W + f^ with 

ii/fiu, <».<-(*>/*>. /?<*^, 

the estimate 

u(0, x n ) , T , / 1 ~ ~ 

sup V " ; <N+-(- sup U + + £ 2 

0<x n <R/2 Xn V -K B r /2,r/2 

remains valid. The proof runs with minor changes. 

Let us compare Theorem 2.3 with results known earlier. Surely, the proof of (ZHO) for 
classical solutions works also for strong solutions if we apply the Aleksandrov maximum princi- 
ple (|AIJ; see also a survey |N05], where the history of this topic is presented). So, it was known 
long ago for b{ G L^. 

In |Li85] the Hopf-Oleinik Lemma was proved for classical solutions of (NDE) in C 1+v 
domains. Note that in this case one can locally rectify dQ using the regularized distance 
( |Li85l Theorem 2.1]). After this Theorem 4.1 |Li85] follows from a particular case b^ = 
of Theorem 2.3, part 2. Similarly, the boundary gradient estimates obtained in [Li86j can be 
reduced to the same particular case of Theorem 2.3, part 1. 

The boundary gradient estimates for solution to (NDE) were established in |LU88] provided 
b G L q , q > n; the Hopf-Oleinik Lemma under the same condition was announced in |NU] . In 
[S10] the second part of Theorem 2.3 is proved for (NDE) under assumption b G L n , b n G L q , 
q > n. In [AN95] the first part of Theorem 2.3 was proved for composite coefficients with 
u(a) = a a , a G]0,1[. 

To compare our result with |S08] , we need an auxiliary statement. 

Lemma 2.4. Let \I/ : [0, <7o] — > be a nondecreasing function. Then there exist nonde- 
creasing C l functions ^ r± : [0, <7o] — > M+ such that < ^ < \& + ; and 

1) if^/(a)/a 2 is summable then (^/ + )'(a)/a is summable; 

2) if^/(a)/a 2 is nonsummable then (\IT~)' '(o~) / 'a is nonsummable. 

Proof. Without loss of generality we can assume <7o = 1. Consider the function ^i(t) = 
\I/(t _1 ) and note that is summable on [1, +oo[ iff ^>(a)/o~ 2 is summable on ]0, 1]. 

Now we define \l/ 2 (t) = L r J )' ( M — r ) + ^i( M )' ( r— L r J )• Using the Cauchy convergence 
criterion, it is easy to check that t\I/ 2 (t) is summable iff is summable. Also it is evident 
that tf 2 (r + 2) < ^(r) < ^ 2 (r -2). 

Finally, we mollify ^ 2 so that ^ 2 (r + 2) < *i(t) < $ 2 (t-2) and set ^ ± (a) = ^ 2 (er _1 =f2), 
expanding to ]|, 1] in a proper way. □ 



In |S08t Theorem 1.8] the Hopf-Oleinik Lemma was proved for solution to (NDE) with 
b = under assumption G dQ and Q D Q^, where 

Q\s = (i G K" : \x'\ < do, <x n < (To}, 

while \l/(cr)/o" 2 is summable. 

By Lemma 2.4, this case can be reduced to Q = Q^+. Then we again rectify the boundary 
and use part 2 of Theorem 2.3. In the same manner, [S08t Theorem 1.9] follows from part 1 of 
Theorem 2.3. 

Remark 4. Note that the assumption (jSj) cannot be removed. Let us describe corresponding 
counterexample (see also [NU] and [S10] ). 

Let u(x) = x n ■ ln a (|x| _1 ) in B^r. Then direct calculation shows that u satisfies an equation 

-Au + b n (x)D n u = with \b n \ < |^n4yzn e L n (B R)R ), 

if R is small enough, and u > = w|a; n =o in £>i?,,,R- However, it is easy to see that D n u(0) = 
for a < and .D n w(0) = +oo for a > 0. 

The condition ()3]) is also sharp. A simple one- dimensional counterexample is given in |A-Z] : 

S t 1 » 

the function 0(s) = J exp ( — / ^p-dr )dt (cf. ([9])) is positive on ]0, 1], vanishes at zero and 
o ^ t T ' 

satisfies the equation 

-4>"( S ) + ^ iff {8) = o. 
s 

However, if u is not Dini continuous at zero then <p'(0) = 0. 

We describe also a more rich family of counterexamples generalizing [S081 Theorem 1.11]. 

Let Q be a convex domain. Suppose that dQ = {x : x n = </>(V)} at a neighborhood of the 
origin, G C 1 , D'<p(0) = 0, and o;(p) = sup |Z)'0(a/)| is not Dini continuous at zero. Let C be 

\x'\<p 

an operator of the form (NDE) in Q with b = 0, and let the condition ([1]) be satisfied. 

It is shown in |ANllj that any solution of Cu = positive in Q and vanishing on dQ at a 
neighborhood of the origin satisfies sup — > as p — > 0. 

Now we rectify dQ at a neighborhood of the origin using the regularized distance and obtain 
a uniformly elliptic operator of the form (NDE) in Br : r with |b(a?)| < 23 "^/-R) f or ^ e 
Hopf-Oleinik lemma fails. 



3 Parabolic case 

In this section we assume d t u, D(Du) G L q ^ t i oc (Q) with some q,£ < oo such that ^ + | = 1. 

We recall an estimate which is a particular case of the statement in \N87\ Sec. 3]. For the 
isotropic case it was proved in |Kr86] . 

Proposition 3.1. Let Ai be an operator of the form (^NDPj in a cylinder Q C BrX ]0,T[, 
and let the condition (TJJ) be satisfied. Suppose that b G L^Q), and a function B such that 
d t M, D(DM) G L 0O (Q) satisfies MM > |b| a.e. in Q. Then for any solution of M.u = f in Q 
satisfying u\ &Q < 0, the following estimate holds: 



The next statement generalizes |AN95| Theorem 2]. For the isotropic case it was proved in 
pOOl Theorem 5.2]. 

Lemma 3.1. Let M. be an operator of the form //NDPj in a cylinder Qr,r, R<1, and let 
the condition (QP be satisfied. Suppose also that the vector function b can be written as follows: 

b = b«+b( 2 ); \b^\eL q4 (Q R , R ), ~ + 7 = l, qJ<oc, (14) 

q I 

\\>^x-,t)\<<& ^ Xn,R) , wGPi. (15) 

Xn 

Then for any solution of M.u = f in Qr : r satisfying w| a ,g < 0, the following estimate 
holds: 

u<N 4 - (|||b«|||^ s ^ + R)< ■ \lf + \U Mu >o h (16) 

where iV 4 depends only on n, v , t, *B and uj. 
Proof. We consider a sequence of operators 

M £ = d t - a ij£ {x; t)D i D j + [b£(x; t) + b%\x; *)]A- 

Here a^ £ are smooth functions satisfying ([T]) uniformly w.r.t. e and tending to a.e. in Qr^r 
as s y 0. Further, 

6^ (a;;*) = max^ 1 ^;*)^" 1 } ■ sign(6, (1) (x; t)); 

hf\ x -t) = [ h ^ X]t ^ Xn>£] 
\bf\x',s;t), x n < e. 

Now we consider the boundary value problem 

M £ v=(^Ul) + l)m + ^^l^-(—Ul)-l) m Qr,r] v\ dlQ =0, 

where p < R will be chosen later. Denote by Be a unique solution of this BVP. By the 
maximum principle ( [Kr76] ). M £ 1} > 0. Define 

i 

253 r 

M £ {x;t) =M i E i) {x;t) + — R / X u (s)ds. 

X n /R 

Then 

M £ M £ = M £ B« + 2<B^ - [&£> + 621 • — ^(^/i?) > |b?>| + |bf| + F(x), 



where 



— 5,(x„/fl)) + as ^^-^ ■ (— 2^(i) - 1) . 

z/ / \ V J + 



x n \ z/ / p 

We set p = min{l; where s~is the root of X u ,(s) = Then, for < a; n < p, we have 



Xr, V V / 



If, otherwise, p < x n < R, then 

F{x) > * «-kM . h _ x„ (1) ) + B ^ . (?» x„ (1) _ A > o. 

o V z/ / o V z/ / + 



So, in any case .M £ B £ > |b £ 1} | + |b £ 2) |. 

Using Proposition 3.1, we obtain the estimate 



u < N(n) ■ ( " j£lkQ + jR ) t ■ \l(M E u) + \l qMu>0} (17) 



v 



for any function u satisfying the assumptions of Lemma. 
Obviously, 

||B 6 ||oo,o< ||B«|| 0OiQ + — RX W {1). (18) 



v 



Further, the function Bg 1 itself satisfies the assumptions of Lemma. Therefore, one can set 
u = Be m f fT7|) and use ( fTHj) arriving at 

Q < 7V(n) • ( lN'Wfi(l + f^(l)) )f . || (A1s B« )+ || aAeiiii (19) 



(we recall that M ( £ 1] > 0). 

If p^HocQ > R{1 + ? 2L(1)) then (USD gives 



^lU.Q + ^l + g X,(l)) /2iV(, v 



{n) 1(1)^ 



(here we used - + | = 1). Thus, in any case we have 

||B«|U,g + R(l + f 21,(1)) „ /2JV(n) .. . ^ „ (lh , V , 2J?(1 + f I„(l)) 



<(^.|| (A4 ,BF) + || a ,«, ea ,„) l + 



z/ \ z/ ' / z/ 

Substituting this estimate into ( II 7p and taking into account the definition of B £ ^ we obtain 
(TTB|) for A^ e instead of M.. Passage to the limit as e — > completes the proof. □ 

The next Lemma is parabolic analog of Lemmas 2.2 and 2.2'. 

Lemma 3.2. Let M. be an operator of the form fNDPj in Q p . p , p < R, and let the 
conditions (T7]) ; [Lfy and < f73]) fre satisfied. Suppose in addition that 

\lb^l\ qAQp , p <m P l (20) 

Let u be a nonnegative solution of M.u = f > in Q P , P - 

1. If u > k on d'Q PtP fl {x n = 0} for some k > £/ien /or £ < | i/ie inequality 



u 



> A: • ^, ^, *B, 21, a;, - ^" 5 (^, ^, ^, ^B, 21, c^, " (p~^ ll|b (1) ||U^, Qp , p + *Ba;(p/i?))) (21) 



holds in Q(i-£Wi-£)p with some positive constants (5 and N 5 . 

2. If u > k on & Q PtP fl {i„ = p} /or some > then for £ < ^ the inequality l[21\) holds in 

\ Q(i-op,cp- 



Proof. We prove the first statement. The proof of the second one is more simple, and we 
omit it. 

First, let £ = |. Consider the barrier function w(x;t) = w(x) + \, where w is defined in 
© with a constant A > ^ n ~ 1 + i_ to be determined later. 

Similarly to Lemma 2.2, direct calculation shows that Aiw < Cg(A, u, Q5)|b^|/7 -1 in Q^-e.. 

Further, w(x; t) < on <9'Q Pi p \ {x n = 0}. Finally, w(x; t) < 1 on <9'Q Pi p n {x n = 0}. This 
gives few — u < on &Q p £_. 

Lemma 3.1, condition (1201) and relation - + I = 1 give for (x; t) <E Q p p. 

u(x;t) > kw{x;t)-kN 4 - (2l* + l)«p? ■ K-M^+H^e^^ > kw(x;t) - kC 10 p-^\lb {1) \\ qAQptP , 

where N4 is the constant from Lemma 3.1 while C\o depends only on n, u, £, A, 03 and 21. 

Similarly to Lemma 2.2, one can choose A = A(n, u, 03, tu) > + 4v ^— j- so large that 

2(1 + A)<p(l/.A) < Then direct calculation gives 

u>k- (^-C r 11 (n,i/,^»,^w>-*|Ib«|, Aepi/J ) + = fe 2 in %_£_x]-p 2 /2,0[. (22) 

Now we consider the set = Q^p , \ Qap P . Note that coefficients b) are bounded on 

H 4 'P 4 >20A 

this set, and 

1 

|||(b( 2 )|||^ p < C 12 {n)Kp-z( j (^P^y ds f < C 13 (n,v,£,<B,u)<Bp?u;(p/R). (23) 

1 

2QA 

We proceed as |LU85| Lemma 3.2] (where the isotropic case was considered) and obtain 

u > k 2 ■ (x(n, v, £, 93, 21, u) - C u (n, u,£, 03,2l,w)p-^|b (1) + b (2) ||| ejt ) in Qp p \ Qp _e_. 

By (T22]) and ([23]) the statement for £ = | follows. 

For arbitrary £ < \ we apply the obtained statement in cylinders Q^^p^p [xP'\ t°) with 
\x°'\ < (1 - 2f)p, (1 - 4£ 2 )p 2 < £° < 0. We arrive at 

u>k- (J^ - C 15 (n, v,£,<8, 21, w) (p~* |||bW |U Qp>p + 95w(p/i2))) + = ^3 

in%_ OA5p x]-(l-3£ 2 )p 2 ,0[. 

Finally, as in the first step, one can proceed as |LU851 Lemma 3.2] in the set <2 PiP \ Q p 
and (EH) follows. □ 

Remark 5. If we replace the assumption / > by /_ e L q i {Q p p ), the estimate (T2"TT) holds 
true with additional term — N 6 (n, z/, £, 03, 21, w, £)p « • ||/-| g /,Q PiP in the right-hand side. The 
proof runs without changes. 

Theorem 3.3. Let Ai be an operator of the form fNDPj in Qr,r, and let the condition 
(Q]) ; (jT^P and [I5\) be satisfied. Suppose also that 

2li(p) = sup P~^\\b m l qA Q p(x o. t0) 0, p^0, (24) 

Qp(x°;t )cQfl,H 



and for p < R 

sup p-^\ AQpAx0 , t0) <^MplR)^ (25) 

Then 

1. Any solution of M.U = f < in Qr,r such that u\ Xn= o < and u(0] 0) = satisfies 

u(0,x n ;0) N+ 
sup < — — • sup u. 

0<x n <R/2 X n R Q r/2 ,r/2 

Consequently, if D n u(0;0) exists then (_D n u) + (0;0) is finite. 

2. Any positive solution of hAu — f > in Qr,r such that u(0; 0) = satisfies 
j n f u(o,x n ,o) > q Consequently, if D n u(0; 0) exists, it is positive. If, in addition, f = 0, 

0<x n <R/2 Xn 

the following estimate holds: 

. nf u(0,x n ,0) > N _ u(0,R/2;-R 2 /2) 

0<x n <R/2 x n ~ ' R/2 

The constants N^ 1 depend on n, v , I, *B ; 93i, 2ti and co. 

Remark 6. If |b (1) | G L q j{Q RiR ) such that q,i < oo and | + | = 1 then ([21]) is obviously 
satisfied. 

Proof. Similarly to Theorem 2.3, we introduce the sequence of cylinders Q Pk ,h k , k > 0, 
where = 2~ k j)g, h k = (kPk, while po < R and the sequence Cfc 4 will be chosen later. 
Denote by M k , > 1, the quantities 

-r?, u(x\t) u(x:t) r~„ . . u(x:t) 

M k = SU P r . > SU P -^ ± - L \ M~ = ^ inf 



Note that in the case 2 M k > 0. 
We define two function sequences 

v k = u- M+h k — ; w fc = M fc /i fe _ - u, 

(p + {h k /R) <p {h k /R) 

where functions ^ are introduced in (jHJ), and denote 

V k = v\\ M k = M+; $ = y? + in the case 1; 
V fc = v\\ M k = M' k : $ = in the case 2. 

It is easy to see that V k \ _ < while the definition of M k gives VJU _, < 0. 

To estimate V k , we consider (x°;t°) G Q Pk ~h k ,h k - Let x° < Then we apply the first part 
of Lemma 3.2 to the function M k h k — V k in Qh k ,h k {x°';t°) (with regard to Remark 5). This 
gives for x G Qh^ ^(x ';^) 

2 ' 2 

Mfc/ifc-FftOr) >M k h k - p(n,z/,£,<B,2li(/i fc ),w,l/2)- 

- AT 5 (n, i/, i, 35, Sti(^), w, 1/2) ■ (Sl^) + »w(flfe/fl))] - 

- N 6 (n, v, £, 03, 2ti(/i fe ) 3 w, 1/2) /i fc * ■ ||(M^ fc )+||, A< 2 hjb , kik (*«>';to). 



By (1211) . we can choose po/R is so small that the quantity in the square brackets is greater 

i 

2 ' 



that ~. As in Theorem 2.3, for (x; t) G Qh^ h^(x°';t°) we arrive at 



V k (x; t) < M k h k ■ 1-13/2 + C 16 (n, u, i, < 3, %i(h h ),u)h. 



i 



In particular, this estimate is valid for (x;t) = (x°;t°). If x° > we get the same estimate 
using the second part of Lemma 3.2. 

Taking supremum w.r.t (x°;t°), we obtain 

sup V k < M k h k • (1 - p/2 + C la mMh k /R)) . 

Qp k -h k ,h k 

Repeating previous arguments provides for m < 

sup V k < M k h k • ((1 - f3/2) m + Cv&i U{ ^j R) ). 

Qp k ~mh k ,h k P/ 2, 

Setting m = L^t^J > we arrive at 

sup Vfc < ^— • exp - A—— + Ci 6 53i — ^- , 

Q Ph+1 , hh 1-/3/2 V V ^7 0/2 J 

where A = — ln(l - 0/2) > 0. 
Therefore, for (x; t) G Q Pk+1 ,h k 



V k (x:t) - — - . 

<M k %, (26) 



max{x n ,/i fc+ i} 

where 7, = ^ ^ ■ ( exp ( - X) + C^B^p). 
Estimate (126]) implies in the cases 1 and 2, respectively, 

where 5^ are defined in ( !T3|) . Similarly to Theorem 2.3, we obtain 

. Tf± k 

Mfc±+1 ^ $7TT ' ± 7i ■ max{l, $(1)}). 

We set Cfc = and choose fco so large and po/R so small that 71 • max{l, $(1)} < |. Note 
that fco and po/-R satisfying all the conditions imposed depend only on n, u, 53, *Bi, 2ti and w. 

Now, as in Theorem 2.3, we observe that the series ^ fc 7fc converges. Therefore, the infinite 
products n =t = Yl(l ± 7^. • max{l, $(1)}) also converge, and we obtain in the cases 1 and 2, 

k 

respectively, 

^ + ^ IP^M* , 1 
fe $(1) 

Thus, all are bounded in the case 1, and all are separated from zero in the case 2. 
Further, we note that M± < ^- sup u. This completes the proof of the statement 1. 

Qr/2,R/2 



If / = then we set K = Qr ; r \ Qr^/2- Similarly to Lemma 3.2, 



sup p t |||b (1) + b {2) i\ qA Q p{x o tt0) 0, p ->• 0. 

Q p (x°,t°)cK 

Therefore, we use the Harnack inequality which can be proved in a similar way as |S10[ Theorem 
3.3] (for bounded lower-order terms see |KrS] ) and obtain 

C 

Aff > -r 11 ■ w(0, R/2; -R 2 /2), 
hi 

where Cyj depends on n, u, £, 93, 2li and to. This completes the proof of the statement 2. □ 
Remark 7. If we replace in the case 1 the assumption / < by / = + with 

sup p-HfPl 9 AQ»,V>W < %Mp/R\ /? < $2 ^l^-, 
Q p Ax '-J KQr,r x n 

the estimate 

sup < iVf • I — sup u + + 

0<x n <R/2 %n Q r/2 ,r/2 

remains valid. The proof runs with minor changes. 

Let us compare Theorem 3.3 with results known earlier. As in elliptic case, the proof of 
Hopf-Oleinik Lemma for classical solutions to parabolic equations with 6j G I ra works also for 
strong solutions by the Aleksandrov-Krylov maximum principle ( |Kr76] ; see also |N05j ). 

In |KHi] the Hopf-Oleinik Lemma was proved for classical solutions of (NDP) in C 1+v, 2 +v 
domains. Using the parabolic regularized distance ( |Li85t Theorem 3.1]) one can locally rectify 
the boundary and reduce the result of [KHij to a particular case = of Theorem 3.3. 

The boundary gradient estimates for solutions to (NDP) were established in |LU88] pro- 
vided Ibl G La+2, Q > n: the Hopf-Oleinik Lemma under condition Ibl G L n r, - + t < 1, 

q,£ < oo, was announced in |NU]. In |AN95] the first part of Theorem 3.3 was proved for 
composite coefficients with b^ G L q+2 , q > n, and u(a) = a a , a g]0, 1[. 

Remark 8. Note that the assumption ( 1251) cannot be removed. Let us describe correspond- 
ing counterexample (see |NUj ). 

Let u(x; t) — x n ■ ln a ((|x| 2 — t)^ 1 ) in Qr,r. Then direct calculation shows that u satisfies an 
equation 

d t u- Au + b n {x;t)D n u = with \b n \ < — — — — G L ^{Qr,r), 

(\x\^ — tp ln((|x| 2 — t)~ L ) 

for any q, £ < oo such that - + j = 1, if R is small enough. By Remark 6, the assumption ( 124|) 
is satisfied. Moreover, u > = m|x„=o in Qr,r- However, it is easy to see that D n u(0;0) = 
for a < and D n u(0; 0) = +oo for a > 0. 

The condition (|T5|) is also sharp. Indeed, considering functions depending only on spatial 
variables we see that the counterexample at the end of Section 2 works also for the parabolic 
operator (NDP). A more rich family of counterexamples also can be extracted from |ANllj . 
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